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Abstract 

The dynamics of massive particles in the collinear high momentum 
regime is investigated. Methods hitherto exploiting large time asymp- 
totic Hamiltonians in the Dirac picture for the treatment of infrared di- 
vergences are adapted to the collinear asymptotic dynamics. The essential 
role of time ordering of the dressing operators is brought out. 
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1 Introduction 



For the treatment of infrared divergencies and the problem of the Coulomb phase it has 
been known for a long time that perturbative contributions up to infinite order have to 
be taken into account. Kinoshita Q, Lee and Nauenberg showed that maintaining 
unitarity in the approximation is the key for solving this problem. Khulish and Faddeev 
@ pointed out that both phenomena are dominated by the asymptotic dynamics of QED. 
They separated out that part of the Hamiltonian which is essential for large time scales. 

The dynamics at large times is dominated by the interaction involving the long wave 
length photons. In the Dirac picture the temporal variation of their contribution to the 
interaction Hamiltonian is small. Actually, for the energy of the photons going to zero the 
corresponding interaction Hamiltonian density in momentum space becomes stationary. 

Similarly, photons collinear to electrons lead to parts in the interaction density sta- 
tionary in the limit of high momentum |Q] . 

In the case of massive particles only in the high momentum limit the stationary regime 
is sustained. 

Adapting the methods of asymptotic dynamics |5|, ^, ^ to the massive case we shall 
show that they lead to Reggeization fs], 10 1. While in the treatment of infrared dynamics 
time ordering of the dressing operator did not come into play, here it will be of essential 
importance. 

We shall also show that the other typical high energy approximation, i. e. eikonalization 
PI , can be obtained using our methods. 

In this paper for developing the appropriate adaptation of our method we study the 
ip^ theory, since it is the simplest model for which Reggeization has been demonstrated. 



2 Approximate Unitarity 

Calculations of bremsstrahlung processes require in a unitary approximation the inclusion 
of an infinity of orders in the coupling constant to ensure a finite total cross section. 

One way of reaching this aim is through an approximate calculation of Dyson's time 
evolution operator [/(t, to) determined by the equation 

i^U{t,to) = Hj{t)U{t,to) (2.1) 
with the well known formal solution 

Uit,to) = Texpi^-ijHi{t')dt'^ , (2.2) 

where Hi{t) is the interaction Hamiltonian in the Dirac picture. 

For definiteness we exemplify our method on the 99^ model. The interaction Hamilto- 
nian 

Hi = -^ f ■.ip^{x):d^x (2.3) 



3! 

can be written as 
9 



Hi = -^ / {v9<-'^(x) + ip^+'>^{x) + 3v9(-)2(x)v9(+'(x) + 3p'--\x)ip^+^^{x)} d^x (2.4) 
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by splitting the field operator ip{x) = + 99'+-' (x) into the creation and annihilation 

parts. 

As in the argument of Kulish and Faddeev Q in the case of bremsstrahlung and 
Coulomb phase we select those contributions to the interaction that are most important 
at large time scales. At first we neglect those terms of Hi that create or annihilate three 
particles and get 

Hi = + H^-'' (2.5) 

with 

= -|y"(^(-)2(a;)<^(+)(a;)d3x , H{-^ = -^J ^^-\x)ip'-+'>^{x)d^x . (2.6) 

or in momentum space 

g f dkdp 



Ht' = -f / -^l^) ""^^P- k)a\k)a{p) exp [,[E{k) + E{p - k) - E{p)]t} (2.7) 



where 



d^k 

dk = - — -.5 — . (2.8) 
(27r)32A;0 ^ ^ 

CoUinear asymptotic dynamics (for p ^ m, where p = \p\) will be incorporated by re- 
stricting the fe-integration to a region ^l{p) in which E{k) + E{p — k) — E{p) tends to zero 
for p ^ 00. K convenient choice for ^{p) is given by the requirement 

kn^ < (2.9) 

where is the momentum vector of the second created particle, = {E[p — k),p — k), 
and cj is a scale parameter. For a > m the region ^{p) defined in this way restricts A;'^ and 
k-^ to an ellipsoid 

'—^ + —<l (2.10) 

with the principal axes a and b given by 

2 p^ (a'^-m?\ a'^-m'^ ,0 cr^ - /„-,-,x 
a = — ^5 ^ H and }f = . (2.11) 

For p ^ a ^ m we have 

a^l . (2.12) 

The calculation of U{t,tQ) is greatly simplified if the time-ordered exponential is suitably 
factorized. This is facilitated by a time dependent similarity transformation. For two 
operators A{t), B(t), not necessary Hermitean, one has 

Texp I (-i) J A{t') dt' + (-i) J B{t') dt' \ = (2.13) 
= Texpj-iy" A(t')dt'|Texp|-i ^ W{^{t' ,to)B{t')Wi(t' ^Iq) dt'^ 
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with 

t 



Wi{t,to) = Te^pl^-i jA{t')dt'j . (2.14) 

The vahdity of the factorization of the T-product follows from a factored ansatz U = WiW 
for the time evolution operator leading to the differential equation 

i^VF(t,to) = W^\t,to)B{t)Wi{tMW{tM (2.15) 
having again as solution a time-ordered exponential 

PF(t,to) = Texp|-iy'H^fi(t',to)5(t')Wi(i',io)dt'| . (2.16) 



For further use we quote a formula already given by Dyson |12] 

W{\t' ,to)B{t')Wi{t\t^) = (2.17) 

OO * *2 

= Y.'"" J "^^nj dtn-l--- J dti[A{ti),[A{t2),...[A{tn),B{t)]...] . 

"=° to to to 

In our context A{t), B{t) will be identified with H^'^^t), H^~\t), respectively. We will use 
the truncated expression 

t 

W{^BWi = B{t) + i j [A{ti),B{t)] dti + 0{g^) . (2.18) 
*o 

Again applying eq. ( |2.13| ) we can introduce 

VF2(t,to) = Texpj -i j B{t')dt'\ (2.19) 



to 

and 

t t2 



VF3(t,to) = Texp|y' dt2 j dh[A{ti),B{t2)\\ (2.20) 
Uo to } 

(note that the time ordering in W3 refers to 12). In this way we get an approximate 
expression for the time evolution operator 

U{t,to) = Wi{tM)W2{tM)W^{t,to) . (2.21) 
It fulfills an approximate unitarity relation of the kind 

?7t(t,to)^7(t,to) =expO(53) . (2.22) 

A similar approach has been successfully used in bremsstrahlung calculations Q. In what 
follows we are going to show that the above approximation produces Reggeization and 
eikonalization of two-body amplitudes for elastic scattering. 
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Let us give the explicit formula for W3 using A{t) = Hi'^\t), B{t) = H[ \t): 

W^{t,to) = Texp^ j d^x j d%(y^-^\yy^^\x)^+{y-x) (2.23) 

with 

A+ix-y) = [ip^+\x),^^-\y)] . (2.24) 

In the expression for W3 we have neglected a part that corresponds to renormalization 
corrections which will not be considered here. 

3 Reggeization 

In the case of soft electromagnetic bremsstrahlung the asymptotic interaction Hamiltonian 
is linear in the photon field. The photon creating part of this Hamiltonian corresponding 
to our H^'^\t) commutes with itself at different times, since the electromagnetic current 
can be well approximated by the classical one. Thus, the time-ordering of the exponential 
is superfluous. 

In contradistinction, H^'^\t), eq. ( |2.6| ), is a product of and (^'+^ and therefore 

does not commute with itself at different times, so that the time-ordering in Wi remains 
essential. A similar comment holds true for 11^2- 

Let us start with a suitable formula for the elastic amplitude Ti{{pi + P2 ^ P3 +P4): 

7if = r + r^=^ , (3.1) 

where 

T=-iJ d^xd'ye-'P^y+'P'^Qix^ - y'){ps,ont\jn{x)jn{y)\piM) (3.2) 

and 

T^' = T{p2^-Va) . (3.3) 



In eq. (|3.2|) the Heisenberg current is denoted by the index "H" . 
The matrix element T can be written as 

T = -iJd'^xd'^ye-'P^y+'P^''e{x°-y^) (3.4) 

x(#3|^7t(xO,-oo)i(x)t/(xO,yO)j(2/)C/(y°,-oo)|pi) . 

In what follows we shall consider small transfers t = (pi — Ps)^ at high energies. It will be 
convenient to work in the laboratory frame of the particle with momentum P2; pi, ps are 
chosen to be large momenta. Detailed calculations will be made in the interaction picture 
and we expect that the leading contribution to the reggeized amplitude is real and comes 
from the quasi-stable multiparticle configurations. Therefore it is reasonable to replace 

U{x',y')^l . (3.5) 
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Inserting a complete set of states between j{x) and j{y) yields 

T = i27r)^5{pi + p2 - P3 - P4)M2 , (3.6) 

with 

M2 = i27rfY.n^n , (3.7) 



n=l 



K = / dfei • . . d^„ ^/f^^^^^,-^-:^/,^^ g:(p3)g„(r) (3.8) 



where 

<^(pi+P2-Er=ife» 
Pi°+P2°-Er=ifc,°+ie 

and 

Gn{p) = {ki,...,kn\j{0)U{0,-oo)\p) . (3.9) 

As we are working in the laboratory frame of the particle with momentum p2 , the current 
j(0) is connected with the interaction of this slow particle. Since we expect that ?7(0, — oo) 
alone accounts for the production of energetic particles, we replace 

.(0)^.(.'->(0).^^'(0) + ^"^(")^^"^(") ) (3.10) 

by 

fiO)=g^^-\0)^^+\0) (3.11) 



in eq. (|3.9| ). 

Next we substitute the approximation for C/(0, — cxd), cf. eq. ( 2.21 ). The operators W2 
and W3 act trivially on the one-particle state \p): 

W2W3\p) = \p) . (3.12) 

So we are left with 

Gn(p) =5(^i,...,4,|(v'<-'(0)v''+'(0))Texp|-i J dtH^^\t)^\p) . (3.13) 

Let us remark that we can limit ourselves to the term of order (n — 1) in the time ordered 
expansion in eq. ( |3.13 ). Lower orders lead to disconnected graphs which should not con- 
tribute to the full amplitude in the asymptotic collinear dynamics considered here. The 
terms higher than (n — 1) necessarily lead to vertex and/or mass corrections neglected 
throughout this paper. 
Therefore 



Gn{p) = (7(-i)""' / dti---dt„_ie(ti,...,t„_i) (3.14) 

—00 

x(0|a(fci) • • • a{K)ip^-\0)ip^+\0)H{''\h) ■ ■ ■ H{+\t^^i)\p) , 
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where the symbol Q{ti, . . . , denotes 

G(ti,...,t„_i) = e(ti -t2)---e(t„_2 . (3.15) 

It is useful to rewrite G„ as 
n 

Gn{p) = ^9:{p) (3.16) 

r=l 

where 



n 
9r 



5(-i)""' / dti • • • dt„„i e(ti, . . . , tn-i) (3.17) 



x{0|a(^i) • • • a{kr-i)a{kr+i) ■ ■ ■ a{kn)v^+\0)H[+\ti) ■ ■ ■ F^+^(t„_i)|p) . 

In evaluating the matrix elements in eq. ( ^.ITf ) we have to commute the a{ki) with Hiitj). 
Thus we shall encounter the following commutators of a{k): 

D{t,k) = -i[a{k),H{+\t)]=ig [ d^x ^p'^-\x)ip''+\x)e''''' , (3.18) 

JxO=t 

d{t,k,k') = [a{k),D{t,k')]=ig f d^xip^+\x)e'^''+''">'' , (3.19) 

[a{k"),d{t,k,k')] = . (3.20) 

The operator d is linear in ip'-^K Therefore we should consider the commutator of d with 
-fTj"''' again. Since d depends on the sum of two intermediate momenta it leads to vertex 
and mass corrections which — as already stated above — are beyond the scope of our work. 
As a consequence, d can be treated as a commuting quantity and it is easy to see that for 
this reason it does not contribute at all. In this way we obtain 

(-i)"-i(0|a(fci) • • • aikr~i)a{kr+i) ■ • • a{kn)y^^+\0)Hi{h) ■ • • Hi{tn-i)\p} = (3.21) 

= Y,{0\^^+\0)D{tl,Kil))---D{tn-l,Kin))\P) , 

where the summation runs over all permutations of the numbers 1,2, ...,r — l,r + l,...,n. 
Applying consecutively 

gi{£;(p-fc)+fcO-pO)t 

D{t,k)\p)=ig \p-k) , (3.22) 

2E[p — k) 

and integrating, we arrive at 

gl^ip) = (3.23) 
= 9 X 

nr 2E{p- k^^ri)) E{p- k^(^n)) + K%) " 

1 1 



2E{p- /c^(„) - A;^(„_i)) E{p- A;^(„) - A;^(„_i)) + /c^^^-, + - P° 



mP- Eh^r K(^) E{p- YX,^r Ki^) + Y.li+r Kii) " P"" 
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Having this form of g", we now look for the kinematical region, that gives a dominant 
contribution. Note that for massive particles the denominators in (7" never vanish. 

We shall calculate the dominant contributions in the leading logarithmic approxima- 
tion, so we consider fixed transverse momenta and 



7r(s) 



» Jkis) + 



» \ Q^J + m' 



for 



s = 1,2,... ,r- l,r + 1, 
and with 



(3.24) 
(3.25) 

(3.26) 



Qr,s =P- ^Ai) ■ 

Then one can show that the terms 
1 



2E{qr,s)[E{qr,s) - 



in eq. (3.23) are different from zero for 



or 



V(s) 



0{ql 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



The second possibility does not contribute in leading logarithmic order, so we are left with 
the conditions eq. ( ^3.29 ). These, together with momentum conservation (see eq. (13.81) ) 



(3.31) 



lead to strong ordering of the A:'s 

pli > A:|(„) » A;[(„_^) » . . . » k\^^^^ » /cl(,_,) » . . . » kl^^^ » A:!, » m . (3.32) 

It means that in the expression eq. ( p.^ ) (with pf = 0) only the same orderings in both 
Gn enter, i. e. 



Vn 



5^" / dA:i---dA:„ 



Pi"+P2"-Er=ife° 
1 



(3.33) 



r = l TTr 



where fl^{'K) is unity in the region defined by eq. ( p.32[ ) and zero otherwise. 
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Integrating over transverse momenta we arrive at 

r=l TVr 



^" = 2fiw {a'Kipr))"' EE4"(vr) (3.34) 



2(27r) 

where 



_ dfcl-..d4 /-(7r)^(pl-E?fcl') 



and 



1 f d^k^ 
K = — — / . 3.36 

2(27r)3 7 (^2 + ^x2)(^2 + (^x_fcx)2) 

The longitudinal integration yields the same result for each I^: 



s (n — 1)! \ \ln s 

where s ~ 2mp\, so that 

5' (ff^Kln^l 

"^-^'W (n-1)! ^ • ^'-''^ 
Then using eq. ( ^^ ) results in 

•^2 = ^E ^ , =gV^+^^ , (3.39) 

n=l ^ ' 

i. e. we have obtained Regge behaviour 

4 Eikonalization 



In this chapter we apply the factorization formula eq. ( p. 21 ) to derive the familiar eikonal 
form of the scattering amplitude. 
We start with 

{p3,Pi\S\pi,P2) , 5" = C/(oo, -00) = VFi(oo, -00)1^2(00, -CX))W3(00, -CX)) .(4.1) 

In this matrix element the part W1W2 can be replaced by unity; other terms of W1W2 either 
give a one-particle s-channel contribution that can be neglected, or vanish identically, 
because W3 conserves the particle number. The first term in the exponent of W3, eq. ( ^.23 ), 
is also non-leading for large c. m. s. energy. 
In this way we arrive at 

{P3,Pa\S\pi,P2) = {p3,PA\Texpi^J r{t)dt + 0{g^)^\pi,p^^ , (4.2) 

where 

t 

r{t) = -g^ j dy^ j d^x d^y tf^-'> {x)ip^-\y)tf^+\x)f^+\y)A+ (x - y) . (4.3) 

-00 x°=t 
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In fact, for large energies the time ordering in eq. ( |4.2| ) is irrelevant. To see this, consider 
any two-particle state \ki,k2) in the center of mass system of ki and k2. Then we have 



r{t)rit')\ku-h) 



-!L J dqdq'a^h + q + q')a\-ki - q- q')\Q) 

{E{ki + q)Y[k^^ - gO - E{ki +q)+ie] 

^ {E{ki + q + q')f[E{ki + q) - q'^ - E{ki + q + q') + ie] 
In the limit ki — > oo, eq. ( [4. 4]) reduces to 



(4.4) 



r{t)r{t')\ki,-ki) ' 



k, >oo 



g 

' A{k 



-oyi j (iq(^q' a){ki+ q + q)a\-ki- q- q)\^) 



q" + q» - le q'^ + g'H - le 

from which it is evident that 

[r(t),r(0]|^i,-fci) ''"="0 . 
As a consequence, eq. ( [1.2| ) can be written as 

r - 

{P3,P'i\S\pi,P2) = {P3,P4\exp { I r{t)dt} \pi,p2) 

with 



(4.5) 



(4.6) 



(4.7) 



r{t) dt \pi,p2) = i 



8(27r)2(pi0)2 7 gi2 + ^2 



d'^q' 



aHpi + q^)a\p2 - q^)\0) . (U 



From eq. (|4.8|) it is straightforward to derive the well-known eikonal formula 



exp^i^o(m|6|) -1 



Mn = 2is J d^e'Pi^ 



(4.9) 



5 Concluding remarks 

We have shown how the collinear asymptotic dynamics described by the Hamiltonian 
eq. ( |2.5| ) restricted to the kinematical region Q, eq. ( |2.9| ), leads to both Reggeization 
and eikonalization of two-particle scattering. While Reggeization is dominated by three- 
particle collinear dynamics, eikonalization corresponds to that kinematical region of the 
asymptotic Hamiltonian where only two particles are collinear. In contradistinction to 
infrared dynamics ^, the time ordering of the dressing operator is essential for three- 
particle collinear dynamics which leads to Reggeization. 
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